Strain-driven isostructural instabilities in epitaxial films 
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Coherent epitaxial growth allows to produce crystalline films with strained structures which are 
unstable in the bulk. Thereby, the relationship between the lattice parameters of the overlayer in 
the interface plane, {a,b), and its minimum-energy out-of-plane lattice parameter, Cmiii(a,fe), need 
not be continuous. This general statement is proven by examples of total energy calculations. As a 
consequence, Cmin, which is determined by the choice of the substrate, and Cmin-dependent intrinsic 
properties of the overlayer cannot always be tuned in a continuous way as one may aim to do by 
means of strained epitaxy. Employing the model of the epitaxial Bain path we predict that such 
discontinuities occur in films of the elements V, Nb, Ru, La, Os, and Ir. The abrupt change of Cmin 
can be exploited to switch properties specific to the overlayer material. This is demonstrated for 
the example of the superconducting critical temperature of a V film which we predict to jump by 
20% at a discontinuity of Cmin. 
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Epitaxy is an important concept employed in the fab- 
rication of film structures with good crystalline quality 
like monoatomic overlayers, multilayers, compound ma- 
terials, and ordered alloys. Such epitaxially grown films 
are technologically important materials with adjustable 
electronic, magnetic, and optical properties. [l| Epitaxial 
growth also allows to study fundamental aspects of low- 
dimensional structures and interfacial effects, and it can 
stabilize structures that arc unstable in the bulk. [2| For 
example, body-centered cubic (BCC) Co does not exist 
as a bulk phase but was stabilized at ambient conditions 
by epitaxial growth on a GaAs substrate. [3] 

In coherent or pseudomorphic film growth, the in-plane 
film lattice parameters are determined by the substrate- 
in-plane lattice parameters. This allows to expose a ma- 
terial to static, non-isotropic strain. Hence, controlled 
homogeneous strain in films is a valuable mean to in- 
fluence their intrinsic properties. [J] This was shown in 
experiments for ferromagnetic bulk- like films, where the 
magnetic moments, the Curie temperature, and the mag- 
netic anisotropy were tuned in a quasi-continuous manner 
by varying the substrate. [5|, Ifil Using piezoelectric sub- 
strates, even a continuous variation of the lattice parame- 
ters and of the related film properties was achieved. [7|, |8[ 

Here, we predict by density- functional calculations 
scanning 30 non-magnetic metallic elements, that epi- 
taxial films of at least six metals show discontinuities of 
their intrinsic properties upon continuous variation of the 
epitaxial strain without changing the lattice symmetry 
(isostructural transition). For example, we find jumps in 
the lattice parameter of La by 0.81 A and in the super- 
conducting critical temperature of V by 20%. 

We model the structure of the considered examples by 
the epitaxial Bain path (EBP), [9[ which describes [10[ 
the bulk part (i.e., the interior) of coherently grown 
body-centered tetragonal (BCT) thick films on substrates 
with four-fold surface symmetry. The EBP is the re- 
lation Cmin (a) between the film lattice parameter in 



the substrate-film-intcrface plane, a = b, and its re- 
laxed (minimum-energy) out-of-planc lattice parameter. 

In the past, the EBP was applied to predict metastable 
states in transition metals, [9|, Illl4l4l | to investigate the 



magnetic order in strained monoatomic overlayers, |15l - 
171 and was used to idei itify the mother phases of strained 
bulk-like films. 0, [H-ig Employing the EBP model, 
we recently predicted three elements, Ru, Os, and U, 
to be ferromagnetic as BCT films in certain ranges of 
strain. ^ 

The examples considered in the present work include 
only the simplest case of non-magnetic, monoatomic 
BCT films. We however emphasize that our main result, 
the possibility of isostructural lattice instabilities 22, l23[ 
under epitaxial strain and a related discontinuous depen- 
dence of intrinsic film properties on the interface lattice 
parameters, is not restricted to certain structures, com- 
positions, or ground states. 

Coherently grown BCT overlayers can be assigned to 
either a BCC, a face-centered cubic (FCC), or a BCT 
unstrained state, which is exposed to in general non- 
vanishing in-plane stress. These unstrained states cor- 
respond to global or local minima of the energy E{a, c) 
in the space of the tetragonal lattice parameters a and 
c. Thus, the unstrained states can be stable phases (cu- 
bic or tetragonal equilibria), metastable states (locally 
stable against arbitrary small distortions), or unstable 
states (only stable against distortions preserving BCT 
symmetry) . 

The absence of forces perpendicular to the substrate- 
film interface allows for a relaxation of c minimizing 
E{a, c) at fixed a. [9|, |lO| Hence, the EBP of a certain 
material is the set of points {(a, Cmin)} C {{a,c)}, for 
which 



EBP 



= I (a, Cmin) E (a, Cmin) = mui £■(«, c)\ . (1) 



The side condition determines the relaxation of c for a 
given a and defines the aforementioned relation Cmin(a). 
Quantities traced along the EBP therefore only depend 

on a, e.g., the total energy E-ebp{o-) = E{a,Caiinio-))- 
The difference between ii'EBp(a) and the ground state 
energy, Eq, corresponds to t he g eneral definition of an 
epitaxial strain energy in Ref. llOl 

We emphasize, that a strict implementation of the no- 
tion 'relaxation into a state of minimum energy', Eq. ([T|), 
leads to relevant, until now disregarded consequences 
concerning heteroepitaxial growth applicable to all ma- 
terial classes to which strained coherent epitaxy applies. 
As we show, the function E{a,c)\a can have two equal- 
valued global minima at c^^^^ ^ c^^^ for certain values of 
a = fldc- Thus, Cniin(a) can be discontinuous (dc) and, 
by epitaxial coherency, also the dependence of Cmin on 
the substrate lattice parameters in the interface plane. 
Provided Cmin(a) is monotonous as usual, it is impossi- 
ble to manufacture films with c e (c^^jj^, c^i^) despite the 
fact that there is no principal restriction on the choice 
of possible substrates. This limits the possibility to tune 
overlayer-material specific properties to a discontinuous 
manner. However, one can switch these properties on 
purpose, if one succeeds to control the lattice parameter 
a in a narrow interval around Odc- 

In the following, we prove the existence of de- 
generate global minima in E{a^c)\a by high-precision 
density-functional theory (DFT) calculations, carried out 
with the full-potential local-orbital scheme FPLO-7.00- 
28. [2J| Using the local-density approximation according 
to Ref. |25| (PW92) and a scalar-relativistic mode for ele- 
ments with atomic number < 49, a full-relativistic mode 
otherwise, we scanned the EBP of 30 elements with the 
atomic numbers 20-23, 29, 30, 38-48^6, 57, 71-80, and 
92 in a wide range of parameters a. |14| The convergence 
of numerical parameters and the basis set were carefully 
checked, [ij] In order to converge the total energy per 
atom at a level smaller than 0.3 meV, linear-tetrahedron 
integrations with Blochl corrections were performed on 
a 24 X 24 X 24 mesh in the full Brillouin zone, apart 
from the elements Cu, Ru, Ir, and Pt, for which a denser 
48 X 48 X 48 mesh was required. 

Vanadium is one of the elements of which the EBP was 
studied several times, [sl, [isl, |26|, |27| This provides us with 
the following, as we will argue incomplete, picture. The 
EBP of V includes two energy minima: one corresponds 
to the stable BCC phase (PW92 equilibrium lattice pa- 
rameter flBcc = 2.929 A); the other one at a = 2.400 
A, Cniin/a = 1.84, is a local minimum, see Fig. [T] This 
local minimum was shown to correspond to an unstable 
BCT bulk state, [l3|, [ij] since its energy is lowered by 
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FIG. 1. Total energy of V with respect to the energy of the 
stable BCC phase, Ebcc = Eq, (left hand ordinate) and the 
DOS at the Fermi level, Ep, (right hand ordinate) for states 
along the EBP. The two discontinuities of the EBP are marked 
by solid bars on both the a and the Cmin/o scale. The first dis- 
continuity at fldc,! ~ 2.637(3) A results in a pronounced drop 
in N{Ep). The second discontinuity at adc,2 ~ 3.345(5) A 
is related to a peak in the electronic DOS. The energies are 
spline interpolated and guide the eye. 



a shear deformation. Tian et al. |19l . |27| however suc- 
ceeded to stabilize coherently grown V-films on Cu{001} 
and Ni{001} which structures could be determined to be 
a strained state of this BCT local minimum state. Ac- 



cording to Refs. % [H, m, and [13, the FCC structure is 
at a saddle point of E{a,c) and a maximum of -Eebp- 
We determined a PW92 lattice parameter of FCC V, 
aFcc = CFcc/V2 = 2.643A. 

It follows from a symmetry argument, that the BCC 
structure {c/a = 1) and the FCC structure {c/a = \/2) 
arc stationary points of E{a, c), [13|, |26| but the side con- 
dition in definition ([1]) determines, whether they are lo- 
cated at the EBP or not. At variance with previous re- 
ports H, [mil [13, we find that the FCC structure does 
not belong to the EBP of V, because it is unstable against 
a tetragonal distortion at fixed a = apcc • This is obvious 
from the plot oi E{apcc,c) in Fig. [2j There, E^apccc) 
has a local maximum at cpcc — v^apcc and the global 
minimum is situated at Cmin = 3.35 A, Cmin/a = 1.27. 

The characteristic feature of this instability is the oc- 
currence of a non-symmetric double well in E(aFcc,c) 
vs. c with the global minimum at Cmin and a further local 
minimum. According to Fig. [2l double wells also occur 
for other values of a in the vicinity of apcc ■ The energies 
of both minima vary with a such that the two minima 
are degenerate at a certain critical value a — a^c- If a 
traverses this value, Cmin jumps in fulfillment of Eq. ([T]) 
by Acmin = |cj„in - c^i^l- — In general, a saddle point in 
E{a, c) (here, at (apcCi v2 apcc)) and a discontinuity in 
the EBP have different a-coordinates, because the double 
wells need not be symmetric. 

In the case of V, the first discontinuity occurs at 
2.635 A < fldci < 2.640 A < apcc- The associated 
change of c^m/o amounts to Acmin/adc,i = 0.2. Wc 
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FIG. 2. (Color online) Cuts of the total energy E{a, c) of V 
for several fixed in-plane lattice parameters a (legend) near 
the FCC saddle point. For each a, the global minimum at 
Cmin is indicated by a filled symbol. Solid lines are spline 
interpolations to the data (symbols) and guide the eye. 



note, that in terms of linear elasticity theory this insta- 
bility of FCC V with respect to a relaxation of c cor- 
responds to a negative elastic constant c^^zz < 0, since 
Czzzz fx d E/dc^. Surprisingly, the EBP of V has yet 
another discontinuity at 3.34 A < adc.2 < 3.35 A with 

ACmin/adc,2 = 0.06. 

Figure [3^ depicts a contour plot of E{a,c/a) of V for 
a large part of the tetragonal parameter space together 
with the EBP. The two minima of E{a, c/a) correspond 
to the two minima of E'ebp from Fig. [T] and the FCC 
structure is the saddle point as demanded by symmetry. 
The contour diagram (see inset) reveals that there are 
pairs of lines of constant energy with E < Epcc which 
enclose one of the two minima and which overlap in their 
a-coordinates in the vicinity of the saddle point. The 
double wells in the plots of E{a,c)\a in Fig. [5] arise from 
this overlap close to the saddle point at apcc- A similar 
peculiarity of the energy landscape at a ^ 3.35 A causes 
the second discontinuity and the related kink of Eebp 
visible in Fig. [TJ 

While the specific energy landscape in the vicinity of 
the saddle point leads to a discontinuous EBP, the mere 
presence of a saddle point — be it symmetry dictated 
or not — does not automatically lead to a discontinuous 
EBP. For example, if abscissa and ordinate axes in Fig.[3K 
were interchanged, determining the EBP would mean to 
search for minima along horizontals, which would result 
in a continuous EBP including the FCC structure. 

In order to analyze the origin of the discontinuities, 
we consider the electronic density of states (DOS) at 
the Fermi level, N{Ep), which is depicted along the 
EBP in Fig. [1] and for the tetragonal parameter space 
in Fig. [3l3. The discontinuity at adc,i is accompanied by 
a pronounced but almost monotonous change in N{E-p) 
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FIG. 3. (Color online) Contour plots in the tetragonal pa- 
rameter space: (a) total energy landscape and (b) N{Ef) of 
V. The two minima of E{a, c/a) are indicated by ■ and the 
saddle point at the FCC structure is indicated by x. The 
inset in (a) shows the saddle point region in more detail. The 
global minimum corresponds to the stable BCC phase, the lo- 
cal minimum is a BCT structure. The EBP (black dot-dashed 
line) has two discontinuities, indicated by black dashed lines. 
Brown hatched lines signal axial ratios of the cubic structures. 
The energy and the DOS were calculated on a dense (a, c) grid 
with a distance of 0.025 A between two adjacent grid points. 



due to a significant change of the unit cell volume. At the 
position of the second discontinuity adc,2j however, there 
is a distinct peak of 1.96states/(eV • atom) in N{Ef) 
along the EBP (Fig. [1]). The same characteristic peak 
appears as tooth-like feature in the large-a region at ap- 
proximately c/a = 0.73 in Fig. [3}d. The EBP of V, which 
crosses this area, becomes discontinuous in order to avoid 
the large N{Ef) (Jahn- Teller behavior). At the local 
maximum of the related double well of i?(adc,2,c) (not 
shown), N{Ef) — 2.4states/(eV • atom). Further anal- 
ysis of the total energy shows that the discontinuity at 
o-dc,i is mainly driven by the kinetic energy contribution. 
Since there is no obvious peculiarity in N{E-p), the dou- 
ble well has to be assigned to an integral effect of the 
whole DOS. 

We searched the EBPs of the 30 elements listed above 
for the occurrence of discontinuous points in the same 



TABLE I. List of elements with discontinuous EBPs in PW92. 
The position of the discontinuity, a^c, the change of Cmin 
across the discontinuity, and the strain energy at Odc are tab- 
ulated. 







discontinuity 


element 


adc 


ACmin 


-EEBp(adc) — Eo 




[A] 


[A] 


[meV/atom] 


V 


2.637(3) 


0.58 


313 


V 


3.345(5) 


0.18 


323 


Nb 


2.933(2) 


0.62 


390 


Ru 


3.010(1) 


0.24 


685 


La 


4.043(1) 


0.81 


178 


Os 


3.063(1) 


0.26 


907 


Ir 


3.076(1) 


0.25 


630 


way as V. 


IJ] The search 


space was 


restricted to 0.7 < 



c/a < 2.1 and had a resolution Aa = 0.05 A. Seven 
discontinuities were found in this way, see Tabic HI Apart 
from V, also the group-5 element Nb (but not Ta) and 
the elements Ru, La, Os, and Ir exhibit discontinuous 
EBPs. For La and Ir, wc identified a negative Czzzz of 
the BCC structure as the driving force of the instability, 
while in Nb the FCC structure has Czzzz < 0. For Ru and 
Os, the discontinuities occur in the vicinity of their BCC 

1 /2 

structures (with cj^i^/adc > !)■ No second instability 
was found for any element except V. 

The present work focuses on non-magnetic elements, 
where we found isostructural transitions under epitaxial 
strain due to either a large peak in N{Ep) or a more 
subtle interplay of all contributions to the total energy 
in the vicinity of a saddle point. In order to find out 
whether the epitaxial relation of a particular substance 
is continuous or not, one usually has to calculate this 
relation in detail. Discontinuities associated with cubic 
Bravais lattices, however, can easily be identified. Since 
these structures are symmetry dictated stationary points 
of E{a, c), [H, [ill they are always on the EBP if Eq. ^ 
is fulfilled. In contrast, an FCC or BCC structure with 
negative Czzzz always results in a discontinuous EBP. 

Similarly, magnetic transitions, like high-spin to low- 
spin I22I or a change of the ground state magnetic or- 
der, [15| arc accompanied by discontinuities of the lattice 
geometry. An example might be the a-dcpendcnce of the 
magnetic ground state of BCT Fe found by Qiu et al. [15| 
These authors did however not discuss the consequence of 
a transition into a different magnetic state for the lattice 
parameter c of coherent films. 

We now return to the case of V to demonstrate, how 
strongly the structural discontinuity may alter another 
particular intrinsic property, namely the superconduct- 
ing critical temperature Tc- Vanadium is a well-known 
strongly coupled superconductor and possesses one of the 
highest Tc-valucs among all elements, T^ = 5.4 K at nor- 



mal pressure. 28| A significant change of Tc may be ex- 
pected at Ode,!, since N{Ep) jumps by 0.2 states/(eV • 
atom) , see Fig. [1] 

The calculations for Tc of BCT V were carried out with 
the ABINIT-DFT package in the version 6.12 [H, HSl 
employing a pseudopotential (PP) generated with the 
FHI98PP code [sij for PW92. With this PP, we found 
slightly shifted values of aecc = 2.963A and of add = 
2.623 A, compared with the FPLO results. The phonon 
properties were obtained through density-functional per- 
turbation theory, [32l | for details about the implemen- 

with a result- 



33 and 34 



tation in ABINIT sec Refs. 
ing electron-phonon coupling constant A = 0.95 and 
a logarithmic-averaged phonon frequency wiog = 288 K 
for BCC V. [35| For comparison, tunneling experiments 
found A = 0.82 ^. We obtained T^ = 3.2 K for BCC 
V according to the Allen-Dynes modified McMillan for- 
mula. [37[ This equation contains one adjustable param- 
eter, the Coulomb-PP /i*, for which we used the conven- 
tional value ^* = 0.13. Across the discontinuity at adc,i, 
A changes by AA = 0.20 and Tc by ATc = 0.6 K or 20%. 
The actual value of Tc may depend on the film thickness 
as it was found in experiments for BCC V films, [38[ but 
we expect that the relative change will not substantially 
be affected by this dependence. 

In conclusion, we have shown that coherently grown 
epitaxial films can exhibit a discontinuous relation 
between in-plane and out-of-plane lattice parameters. 
Thus, it is impossible to manufacture films with an out- 
of-plane lattice parameter which falls into the unstable 
region, although there is no principle limitation on the 
choice of the substrate. In particular, we predict that 
BCT bulk-like films of the elements V, Nb, Ru, La, Os, 
and Ir undergo isostructural lattice transformations in- 
duced by epitaxial strain. These lattice transformations 
are necessarily accompanied by abrupt changes of all in- 
trinsic properties of the films which continuously depend 
on the structural parameters. As a consequence, the film 
properties can be switched, if one succeeds to control the 
substrate lattice parameter appropriately. We demon- 
strated this possibility by the example of V, where the 
superconducting critical temperature changes by about 
20% at a critical strain value. 

Since the mechanisms of the discussed instabilities may 
occur in arbitrary systems, our results are not limited 
to monoatomic systems, BCT films or specific ground 
states. Verification of the predicted lattice transforma- 
tions should be possible by advanced diffraction methods 
applied to coherently grown thick overlayers on a series 
of substrates with quasi-continuous lattice spacings. 

Wc thank K. Doerr, S. Gemming, S. Haindl, S. Fahler, 
and A. Mobius for discussions. 
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